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ABSTRACT 

We show that dwarf nova disks in quiescence have rather low magnetic Reynolds number, of order 
10 3 . Numerical simulations of magnetized accretion disks suggest that under these conditions magne- 
tohydrodynamic turbulence and the associated angular momentum transport is sharply reduced. This 
could be the physical origin of episodic accretion in dwarf nova disks. If so, the standard disk instability 
model needs to be revised. 

Subject headings: accretion, accretion disks — MHD — turbulence — novae, cataclysmic variables 



1. INTRODUCTION 

It is a commonplace that resistive diffusion of magnetic 
fields ia negligible in aatrophyaical plaamaa. — The impor 



tanco of resistive diffusion can be measured via the mag 



netic Reynolds number 



Re 



M 



LV 



where rj is the resistivity, and L and V are a typical 
lengthscale and velocity, respectively. Reu is often large 
in astrophysics because the length and velocity scales are 
large. For example, the decay time of magnetic fields in a 



sunspot, as estimated in the introduction to Parker (1979) 
is 300 yrs, implying a magnetic Reynolds number of order 
10 7 . The decay time of magnetic fields in the center of the 
sun, or in the galactic disk, is even longer. In this letter we 
discuss an uncommon situation where Reu is small and 
the decay time is short: dwarf nova disks in quiescence. 

Dwarf novae (DN) are binary stellar systems in which 
a white dwarf accretes matter through an accretion disk 
from a main-sequence companion. DN go through regular, 
but not periodic, outbursts in which the luminosity of the 
disk increases sharply. They exhibit a rich and com plex 
phenomenology, reviewed in detail by Warner (1995). DN 
are often regarded as "laboratories" for the study of as- 
trophysical disk systems because they are bright, nearby, 
and have short characteristic timescales; changes in our 
understanding of these systems therefore have broad im- 
plications for other disk systems such as X-ray binaries, 
young stellar objects, and active galactic nuclei. 

The significance of low ReM in DN disks is connected 
to recent advances in the theory of angular momentum 
transport in disks. Balbus & Hawley (1991) showed that 
disks are linearly unstable in MHD, thereby revealing the 
long-sought-for instigator of disk turbulence. More recent 
work on convection and nonlinear hydrodynamic instabil- 
ity makes it seem unlikely that these once promising al- 
ternatives play any role in angular momentum transport 
in disks. Simulations of the nonlinear development of the 
Balbus-Hawley instability show that it transports angular 
momentum outward and acts as a dynamo in the sense 



that it su stains a magnetic field in the presence of dissipa- 
tion (see Balbus & Hawley (1997) fo r a complete review). 
Most signific antly, simulations by Hawley, Gammie, & 



Balbus (1996) (hereafter HGB) show that the nonlinear 



development of the instability is sensitive to the presence 
of resistive diffusion. A series of experiments in that paper 
shows that at Reu = L 2 z ^l/rj — 10 4 the magnetic field is 
depressed, 2 while at Reu = 2000 the magnetic field dies 
away, as does the angular momentum flux. This suggests 
that at Reu <, 10 4 MHD turbulence and the associated 
angular momentum transport may die away. In this letter 
we show, using SS Cygni as a specific example, that Reu 
is of order 10 3 in dwarf nova disks in quiescence. 

In §2 we briefly summarize the standard disk instability 
model for DN outbursts. In §3 we estimate ReM for SS 
Cygni's disk in quiescence. In §4 we tentatively put for- 
ward a new limit cycle model for DN outbursts. We discuss 
the implications of these results in §5 and conclude in §6. 

2. THE DISK INSTABILITY MODEL 

Historically, a variety of mechanisms have been pro- 
posed to explain DN. The disk instability model, invoking 
an intrinsic modulation of the accretion rate in the disk, is 
now generally accepted as the explanation of DN outbursts 
(see Cannizzo (1993b)| for a review and a discussion of the 
historical development of the subject). The essence of the 
mechanism is that the disk cycles between two states, a 
large accretion rate state (hot and mostly ionized) and 
a low accretion rate state (cold and mostly neutral). The 
disk cannot settle into a steady intermediate state because 
that state is thermally unstable. 

The simplest version of the disk instability model, in 
which the a ngular momentum transport efficiency ct (in- 
troduced by [ghakura fc Sunyaev (1973) ) is held co nstant 
throug hout the cycle, does not work, however (e.g 



Smak 



(1984)). Detailed time-dependent calculations show that, 



in order to obtain an outburst of sufficient amplitude, a 
must vary between the hot and cold state. Researchers 
have been able to reproduce outburst light curves and 
mean times between outbursts using ahot — 0.1 and 
a co id — 0.01. Another oft discussed possibility is that 
a obeys a scaling relation such as a = cto(H/r) q , where 



1 also DARC-CNRS, Observatoire de Paris-Meudon, 92195 Meudon, Cedex, France. 
2 Notice that a different definition of ReM is used in HGB. 
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H is the disk scale height and r the local radius. To our 
knowledge this prescription has been used in the DN con- 
text only to study decay properties of the outburst light 
curves. Both these prescriptions for varying a are phe- 
nomcnologically motivated; while there exist some post hoc 
theoretical justifications for varying a, in our view they are 
not strongly physically motivated. 

3. RESISTIVITY IN SS CYGNI 

What then is Rcm in dwarf nova disks? For dcfmitcncss, 
we focus on a nearby, well-studied system: SS Cygni. SS 
Cyg has a 1.2 M Q white dwarf accreting from a disk fed 
by a 0.7 M K5V companion, a period of 6.6hr, a distance 
from the primary to the LI point of about 6 x 10 10 cm, a 
mean interval between outbursts of 40d, and a characteris- 
tic decay time from outburst of 2.4d. The mean accr etion 
rate i s not we ll known but is of order 10 -9 M Q yr _1 ( Can 
(1993a)| ) 



We require a temperature and density to estimate the re- 
sistivity in SS Cygni in quiescence. These are only weakly 
constrained by observations. In the absence of direct mea- 
surements, we must turn to a theoretical model for the 
disk evolution: the standard disk instability model. Fu- 
ture observations may provide better constraints on phys- 
ical conditions in the disk. 

SS Cyg has been theoretically studied in detail using the 
standard thermal limit cycle model (see Cannizzo (1993a) 
for the most detailed study to date). We have run our 
own evolutionary models of SS Cyg to estimate physical 
conditions in the quiescent disk, using a modern, time- 
dependent, implicit, adaptive grid code flHamcury et a.lj 



1997 



We adopt the same parameters for SS Cyg as 



Cannizzo's standard model (ahot = 0.1, a co id — 0.02), and 
our calculated £ and T c are similar to those in his stan- 
dard model. At a radius of 2 x 10 10 cm from the primary 
a typical surface density in quiescence is 200 gem -2 , and 
a typical central temperature is 3000 K. Then £1 = 4.5 x 
10" 3 s- 1 , c s = 3.5 x 10 5 cms-\ H = c s /Vl = 7.7 x 10 7 cm, 
and p ~ £/(2i2") = 1.3 x 10 -6 gcm -3 . In LTE hydrogen 
is then predominantly molecular. The disk is marginally 
optically thick in a Rosseland mean sense. 
The resistivity is given by 



V 



c 2 rn e v c 



where v c is an effective collision frequency for electrons. 
Since electron-neutral collisions dominate, we take v c = 
n n o~ en v, where n n is the neutral number density, o~ en is the 
electron-neutral momentum exchange cross section, and 
v = x /l28fcT/(9 7rm e ). At 3000 K er e „ = 1.3 x 10 -15 cm 2 
flHayashi (1981)1 ). Then rj = 1.67x 10 4 n„/n e . The electron 
abundance at this temperature is determined primarily by 
ionization of Na, and to a lesser extent, Ca and K (we 
assume solar abundances). A detailed solution for ion- 
izatio n equilibrium (we have used a code kindly provided 



to us |by Dr. P. Hoeflich, but a simplified calculation in- 
volving Na, C, and K gives nearly identical results) gives 



3.73 x 10 17 cm~ 3 . and n e 



j.53 x 10 n cm~ 3 . so 



7 29 x 10 9 cm 2 s- 1 . 



V 

The natural length and velocity scales in disks are H 
and c s , respectively. The disk magnetic Reynolds number 
is then Rcm = c s H/r). Using the resistivity just calcu- 



lated, Rcm = 3670. This is below the value at which 
MHD turbulence is depressed in the simulations of HGB. 
Similarly low values of ReM are obtained at other radii 
in our model. Of course, this Rcm was obtained assum- 
ing a co id — 0.02. At lower a the disk will be cooler and 
denser, and then ReM will be smaller since r\ is very sensi- 
tive to temperature. Since a depends on the resistivity as 
well, there could be a runaway decay of MHD turbulence. 
All this suggests that, at least in SS Cygni, angular mo- 
mentum transport may die away once the disk goes into 
quiescence. 

Before going on, let us briefly consider two other issues. 
First, does ambipolar diffusion play any role? The rela- 
tive importance of ambipolar diffusion and resistivity is 
controlled by the product of the electron and ion Hall pa- 
rameters, DR = Lo e LOij (venVin), where u) e j are the ion and 
electron cyclotron frequencies and Vi n is the ion-neutral 
collision frequency; when DR > 1 ambipolar diffusion 
dominates. Assuming equipartition magnetic fields at our 
fiducial point in the SS Cyg disk, we find DR = 0.02. 
Since the field is likely to be weaker than this, resistivity 
dominates ambipolar diffusion. 

Second, can nonthermal ionization save the day? Here 
the chief concern is X-rays; dwarf novae typically have 
L x ~ 10 30 ergs _1 , but SS Cyg is particularly bright in 
hard X-rays in quiescence, with Lx (l — 37keV) ~ 1.5 x 
10 32 ergs _1 ( |Yoshida et al. (1992)| ). Only photons with 
E ;> lOkeV will penetrate the disk. These will Comp- 
ton scatter and diffuse downward into a layer of thick- 
ness ~ 30g cm~ 2 corresponding to an effective optical 
depth of 1 (Glassgold et al. (1997), [gea & Glassgold 
(1997)). The intercepted flux per unit area of the disk 
is f x ~ L x (H/r)/(4irr 2 ). Suppose this flux diffuses into a 
layer of thickness H and produces 1 ion per Ei = 36 eV. 
Then the volume ionization rate is £ = L x j (Aivr^Ei) cm -3 . 
Balancing this against dissociative recombination at a rate 
n\ x 8.7 x 10" 6 T- 1 / 2 cm~ 3 , we find n e = 1.1 x 10 9 cm~ 3 . 
Since this is much less than the LTE electron number den- 
sity we can neglect X-ray ionization in the bulk of the disk. 
It is possible, however, that a thin layer on the surface of 
the disk will be ionized by X-rays, and accretion will pro- 
ceed in that laye r in a manner similar to that described by 



Gammie (1996) in the context of protoplanetary disks. A 



thin partially ionized l ayer could explain eclipse maps o f 
quiescent dwarf novae (Wood et al. (1992), Home (1993)), 
which are interpreted as showing emission from a warm 
optically thin disk. 

Suppose that MHD turbulence dies away in SS Cyg's 
disk in quiescence. It is possible that some weaker resid- 
ual transport process will be present in the disk. Let 
us consider some of the possibilities. Convective turbu- 
lence is ruled out for at least two reasons: the disk is 
only marginally optically thick, so convection will not be 
present, and even it were present, recent studies show that 
convective turb ulence produces inward angular momen- 
tum transport ([Ryu fc Goodman (1992) , [Stone & Balbui 



(1996), Cabot (1996)| ). Turbulence due to nonlinear hydro 



dynamic instability also see ms unlikely in light of recent 
analytic and num erical work (Balbus et al. (1996), Hawley 



& Balbus (1997)) which shows that nonlinear instability 
is simply not present in a Keplerian disk. Gravitational 
instability is ruled out, since for reasonable values of the 
surface density and of the "floor" temperature, set by irra- 
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diation of the disk by the secondary, hot spot, and primary, 
the disk has Q ^> 1. The remaining possibilities involve 
a spiral wa ve or shock, driven either by the comp anion's 



frequency. The results are shown in Figure 2. A fit accu- 
rate to 5% is 



a spiral wa ve or snock, driven eitner by tnc comp anion s 

tidal field ( [Sawada et al. (1986)1 iSpruit (1987)1), or by log(T 4 ) = 3.867 + 0.038 log(S) + 0.20 log(fi) + 0.01 log(fi) 5 



a global linear instabil ity (Papaloizou & Pring h i [mi 



Goldrcich et al. (1986)| ) 



The role of tidally driven spiral shocks is controve rsial 
(Rozyczka & Spruit (1993), (gavonije et al. (1994) ; see 
also the cla imed detection of spi ral waves in IP Peg in 
outburst by Steeghs et al. (1997) ). In our view, Savonije 
et al. have shown that a disk with small H/r and fixed 
mass couples only weakly to the tidal potential, and suffers 
no global instability. They have not, however, considered 
the interaction of the disk with the mass transfer stream. 
Since the incoming material has approximately constant 
specific angular momentum, which will be conserved in 
the absence of turbulent diffusion, it will accumulate in a 
partially press ure supported ring. This ring or torus may 
then suffer the [Papaloizou fc Pringlc (1984) i nstability, or 
an al lied glob al instability, as suggested by Rozyczka & 
Spruf j (1993)1 

4. A NEW LIMIT CYCLE? 

The possibility of global hydrodynamic instabilities in 
DN disks leads us to tentatively propose a variant of the 
classical limit-cycle model for DN outbursts. In the hot 
state the evolution proceeds as usual until accretion causes 
the surface density to drop below the point where the disk 
can maintain the hot solution. The disk then cools, re- 
combines, and forms molecular hydrogen. If the magnetic 
Reynolds number is low enough, then MHD turbulence 
will decay away. The disk then enters quiescence (the vis- 
cosity could even be as low as the molecular value, if no 
other angular momentum transport process is at work), 
and new material accumulates in a partially pressure- 
supported ring in the outer disk. Eventually this ring 
suffers a global hydrodynamic instability. The resulting 
shocks raise the central temperature of the disk, and hence 
ReM: above a critical value. MHD turbulence sets in, rais- 
ing the temperature still further and providing significant 
angular momentum transport. An outburst results. This 
idea is illustrated in Figure 1. 

It remains to be seen whether this scenario can be in- 
tegrated into a successful time-dependent model for the 
outburst. We suspect that at least one additional ingre- 
dient is required, and that is a modification of the classi- 
cal di ffusion equatio n for surface density in a viscous disk 
(e.g. Pringlc ( 1981 )| ) . An important physical assumption 
underlying this equation is that the local stress adjusts in- 
stantaneously to a value specified by the current s urface 



where £ and fl are expressed in gem -2 and s _1 respec- 
tively. Notice that, for a given value of a and for realistic 
opacities not all regions of this plot are accessible. In par- 
ticular, disks in the upper right corner are unlikely to be 
cool enough to reach T 4 . 

The above discussion is predicated on the idea that 
MHD turbulence in disks decays below a critical ReM of 
order 10 4 . This is consistent with the evidence from nu- 
merical simulations, but cannot be regarded as firmly es- 
tablished given the small number of simulations that have 
been done to date and their low numerical resolution. 3 
One possible complication is that saturation may depend 
on magnetic Prandtl number P tm = vjr\ rather than sim- 
ply on Rcm (Batchclor (1950)); here v is the microscopic 
viscosity. The sense of the theoretical argument is that 
when viscosity is larger than resistivity the field builds up, 
while in the opposite limit the field dies away. Counterin- 
tuitivcly, this argument suggests that a should increase as 
v increases. There is weak evidence for this from the simu- 
lations in that those with larger artificial viscosity saturate 
with slightly larger a. Since microscopic viscosity in disks 
is small compared to the numerical and artificial viscosi- 
ties present in the simulations, there is the possibility that 
MHD turbulence may decay at even higher temperatures 
than we have indicated. Clearly there are many unan- 
swered questions about the saturation of disk turbulence 
that might profitably be addressed by future numerical 
simulations, with the promise of direct application to as- 
tronomically interesting systems such as DN. 

Several characteristics of the new limit cycle are worth 
mentionin g here. First, our model b ears some resemblance 
to that of Armitage et al. (1996), who speculated that 



MHD turbulence would cease as DN enter quiescence be- 
cause the magnetic field would become superthermal. Our 
physical mechanism is, however, different and less specu- 
lative. We expect the field to resistively decay, rather than 
being frozen in and dynamically dominant. Second, in our 
scenario the outburst is triggered rather differently than 
in the disk instability model. Assuming that turbulence 
exists above a critical ReM outburst should begin when 
the disk rises above a critical temperature, not when the 
disk becomes thermally unstable. Third, the details of the 
outburst will be sensitive to metallicity, or at least the 
abundance of Na, because Na + provides most of the free 
electrons at low temperatures. Finally, preliminary calcu- 
lations suggest that the model is applicable to X-ray bi- 
naries, although disk irradiation could complicate matters 
somewhat. 

6. CONCLUSION 

We have argued that MHD turbulence and the associ- 
ated angular momentum transport may die away in dwarf 
novae in quiescence. If there are no other significant 
sources of angular momentum transport in disks, then this 
is the physical origin of episodic accretion in dwarf novae. 

3 The critical value of the resistivity, if any, must be determined from nonlinear theory. Linear theory says that stability is recovered at a 
value of the resistivity that depends on the field strength. Absent an almost thermal external field, however, the field strength in the disk is 
determined by the nonlinear outcome of the instability itself. 



density and temperatu re. Simulations of HGB and Bran- 



denburg et al. (1995) suggest that there is a relaxation 



time of many orbital periods as the stress adjusts to the 
current equilibrium value. To our knowledge, no study of 
time-dependent behavior in disks includes this effect. 

5. DISCUSSION 

As a guide to where resistive effects become important 
in disks, we have calculated the temperature T4, where 
Rgm = 10 4 , as a function of surface density and rotation 
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We have proposed a scenario in which a global hydrody- 
namic instability heats the outer disk, thereby raising the 
conductivity and initiating the outburst. This scenario is 
a modification of the standard disk instability model that 
provides a physical explanation for episodic accretion yet 
retains many of the standard model's most attractive fea- 
tures. 
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Fig. 1. — A schematic of the proposed outburst cycle. The solid line ("S-curve") shows the thermal equilibrium solution for the disk in 
SS Cygni at radius r = 2 X 10 10 cm, using the standard assumption of a = 0.1 on the hot branch and a = 0.02 on the cold branch (the 
transition temperature is 2.5 X 10 4 K). Our estimates show that when the disk "falls off" the hot branch, it will pass through the line marked 
ReM = 10 4 , where MHD turbulence dies away. The disk is then passive, and is reactivated only by a global hydrodynamic instability. 



6 



Episodic Accretion in Dwarf Novae 




10" 7 10" 6 10- 5 0.0001 0.001 0.01 0.1 1 

Qts" 1 ] 

Fig. 2. — The temperature T4 where ReM = 10 4 , as a function of rotation frequency and local surface density. The heavy solid lines are 
spaced at intervals of 10 3 K, the light lines at intervals of 200 K. Recall that at T > T4, Rcm > 10 4 . 



